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Abstract 

The radiation reaction fields are calculated for an accelerated changing dipole 
in scalar and electromagnetic radiation fields. The acceleration reaction is 
shown to alter the damping of a time varying dipole in the EM case, but not 
the scalar case. In the EM case, the dipole radiation reactionfield can exert 
a force on an accelerated monopole charge associated with the accelerated 
dipole. The radiation reaction of an accelerated charge does not exert a torque 
on an accelerated magnetic dipole, but an accelerated dipole does exert a force 
on the charge. The technique used is that originally developed by Penrose for 
non-singular fields and extended by the author for an accelerated monopole 
charge. 

It is well known that an accelerated charge radiates and that the emitted field from the 
accelerated charge exerts a force ( the radiation reaction force) back on the charge itself. 
This radiation reaction force is usually derived either by an appeal to a balance of the 
energy and momentum emitted by the charge or by a detailed examination of the energy- 
momentum tensor just near the charge. These techniques invariably have difficulties with the 
fact that the fields from a point charge diverge, and necessitate re-normalisation of various 
quantities (eg the mass) in order to extract reasonable results. However, it was discovered 
by Penrose and Unruh that the radiation reaction field for a charge could be extracted from 
the radiation field by means of an integral over the future null cone of the the particle. This 
integral, which in the absence of any sources, exactly gives the field strength at a point, also 
gives a finite result whan applied to the field emanating from a point charge, and that finite 
result is exactly the radiation reaction fields at the particle. In this paper I will show that 
this approach can be generalised to the case of an accelerating dipole and leads to finite 
radiation reaction fields at the location of the point dipole. These fields lead to the damping 
forces on that accelerating dipole. These results will ultimately be applied in another paper 
on the question of the equilibrium polarisation of an accelerating particle with spin, but 
seem to be of sufficient general interest that they are here separated out for detailed study. 

In this paper, I will derive the radiation reaction fields at the location of point source for 
a massless scalar field, and the reaction fields at the location of a point dipole source for the 
electromagnetic field, using the Penrose integral to do so. Thus I will begin by giving a very 
brief review of the Newmann-Penrose (NP) spinor formalism, and the natural null tetrad 
metric for an accelerated path in Minkowski spacetime. That accelerated path is assumed 



for the purposes of this paper to consist of an acceleration restricted to a plane (linear or 
circular acceleration, but with non-constant acceleration). 

I. NULL METRIC, TETRAD AND NP SPIN FORMALISM 

This section will be a very quick review of the NP spinor formalism, presented as much 
to specify the notation I will use as for any other purpose. For a more complete introduction 
to spinors, see Penrose and Rindler [0] 

Given a path in spacetime, define a ur8<j) coordinate system as follows: 

Parameterize the curve by its length parameter u. (ie, ^/^^f - ^~ = ^ each point 
along the curve, define the future directed null cone centered at that point along the curve, 
and label that null cone by the parameter u. 

On each of the null cones, choose the direction of the acceleration vector, and call it the 
9 = direction. Along each of the null generators of the cone (null lines which originate 
at the point along the curve at the vertex of the cone) define a radial affine parameter r, 
such that the tangent vector to the null curve parameterized by r (which is a null vector) 
has a dot product with the tangent vector to the original curve of unity- Ie, n^^r 1 = 1. 
Now, the two spheres defined by u and r constant are metric two spheres on which we 
will define angular coordinates 9 and 0. They are chosen so that the metric on these two 
spheres is the usual two sphere metrics.. Because of the way that r has been defined, the 
circumference of these two spheres will be r so the metric on these two spheres will be 
r 2 (d8 2 + dip 2 ). The direction 9 = is at the null generator pointing in the direction of the 
acceleration vector, and 9 and <fi will label the null generators. This procedure will define a 
new coordinate system centered on the path of the particle. The metric of flat spacetime in 
these coordinates is 

ds 2 = (l + 2f(u)r cos(#) - r 2 (g(u) cos(0) + f(u) sin(fl)) 2 - r 2 g{uf sin(0) 2 cos(#) 2 ) du 2 

+2dudr + 2{f{u) sin(6 ) ) + g{u) cos(0) )r 2 dud9 — 2g(u) sm(9) cos(9) sm((j))r 2 dud(f) 
-r 2 d9 2 -r 2 sm(9) 2 d(f) 2 (1) 

where f{u) and g{u) will be related to the acceleration. f{u) is the acceleration, and g{u) 
is the rate of change of the direction of the acceleration. The acceleration is assumed to be 
confined to the place <fi = or tt. 

We now define a null complex tetrad for this metric. This is a set of four vectors Z M , 
m M and # such that each of these vectors is null (has zero inner product with itself) and 
such that 

Pn M = m^m^ = 1 (2) 

and all other inner products are zero. The vector n M has already been defined as the vector 
tangent to the null generators of the null cones. The vector m M is defined as a complex 
vector lying tangent to the two surface of the u and r constant spheres. m M is the complex 
conjugate to fi. Thus the null vector Z M will be orthogonal to the these two spheres (as is 
n M ). Z M is assumed to be such that is greater than zero (ie, is a future pointing null 

vector.) For definiteness, I will choose = + ie'i where eg and are the unit vectors 
which lie along the u, r, <ft constant line and u, r, 9 constant lines respectively. 
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In the above ur9<j) coordinate system, these vectors therefor are 



n"= [0,1,0,0] (3) 
1 



y/2 
1 

75 



0,0,-, 



r ' r sin(0) 

1 -i 
0,0,-, 



r ' r sin(0) 



(4) 
(5) 



p=[l,-I-/(u)rcos(0), /(u)sin(0)+^(«)cos(0), - </(u) cot(0) sin(0)] (6) 

In addition to the null vectors, the formalism defines a set of two dimensional complex 
spinors. In particular they assumed that there are two separate two dimensional spinor 
spaces which are anti-unitarily related to each other (these are just the two unitarily in- 
equivalent spin 1/2 representations of the Lorentz group). Tensors over these two spaces 
are designated by indexed symbols whose indices are capital Roman letters. Tensors over 
the one representation have plain indices and tensors over the other are designated with 
primes on their indices. Because of their anti-unitary relationship, there exists a mapping 
from one type of tensor to the other, which we will denote by complex conjugation. Thus 
(w A )* = w A . This is defined so that the inner product of complex conjugate vectors is 
just the ordinary complex conjugation. (w a Va)* = w a Va'- A complete set of basis vec- 
tors on these spinor spaces are designated by l a , o a in the one case, and i A \ o A ' in the 
other case. They are chosen so that (i A )* = i A \ etc. These spinor spaces are related to the 
spacetime vectors by means of the spin matrices cr^ A ,, matrices which for a given vector 
represent a mapping from the one spinor space to the other via V^cr^A'- H we choose the 
basis vectors appropriately in the spinor space, then these matrices are just the four Pauli 
spin matrices, 1, cr x ,a y ,a z . In particular I will assume that they are chosen so that 

l^ = (l + a z )/2 (7) 

= (1 - a z )/2 (8) 

= (a x + i * a y )/2 (9) 

m M cr M = a x - i * <r y )/2 (10) 

with l a and la' both represented by (0, 1) and o^ and oa 1 both represented by (1, 0). Thus, 
we have 

^°A4'= °AOA> (11) 

n^CTAA'= l^A^A' (12) 

m ^ (J AA'= °AtA> (13) 

"V°A4'= laOa' (14) 

From now on I will freely alternate between the spinor representation and the vector 
representation. Thus a tensor S^ u can also be written as Saa'bb 1 = S^uO'aa'O'bb'i e ^ c - 

These spinor spaces also have a metric defined on them, a metric which must be compat- 
ible with the Lorentzian metric. The metric, designated by either e AB for the spinor, or e A B 
for the anti-spinor space is an antisymmetric metric, such that e AB = —e BA . The inverse 
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metric is eab which is also antisymmetric and is chosen so that €ab£ CB = Indices are 
raised and lowered by means of the metric, but because of the antisymmetry of the metric, 
the order of the indices is of crucial importance. My convention, following Penrose and 
Rindler §, is 

^ = e AB uj B = -u B e BA (15) 
u A = uj B €ba- (16) 

and similarly for the prime space. The basis vectors i A and o A are chosen to obey oal a = 1, 
which is compatible with the relation of these spinor bases to the spacetime null tetrad, and 
the normalisation of those tetrad vectors. 

This antisymmetry of the metric and the attendant risk of confusion in index manip- 
ulations (together with the presence of the two separate types of spinors) is probably the 
greatest impediment to the adoption of spinors as a standard approach in special and general 
relativity. Given these disadvantages, there must be some advantages to the spinor formal- 
ism which would persuade any but masochists to adopt the formalism. The key advantages 
are firstly that spinors allows a unified treatment of the various fields of whatever spin type 
in one simple notational system. Ie, scalar, spin 1/2, vector, spin 3/2, .... fields can all be 
treated very similarly, in an extremely compact and transparent way. The second advantage 
is that spinor space is a two dimensional space. This means that there is only antisymmetric 
tensor of rank two, and all other antisymmetric tensors of rank two must be proportional to 
this tensor. Furthermore, since the metric is antisymmetric, it can be chosen as the fiducial 
antisymmetric tensor. Any two indices (of the same type) of a spinor tensor can always be 
written as a combination of symmetric and antisymmetric pairs. Thus, we can write any 
tensor 

S ...A..B.. = S ...(A\..\B).. + 1 CDS ....a.D.. e AB r IT) 

Any tensor can thus be written as the sum of tensors which are totally symmetric on their 
indices times products of the metric tensors. This ability to represent any tensor as either 
entirely symmetric tensors or multiples of the metric is the key power of the spinor notation, 
and achieves its greatest power in the representation of massless fields. A massless field of 
spin s is represented by a tensor with s indices, all of which are completely symmetric, say 
^ab..s where the tensor is symmetric under interchange of any two indices. Furthermore, 
the equations of motion of a spin s massless field are simply written as V AA ^ ab...s = 0, 
where V A A< is the covariant derivative, defined on spinors such that a\ A , and the metric tAB 
are covariantly constant. It is the compactness of the spinor notation, and the transparency 
of the symmetries of the fundamental tensors which give the spinor notation its power. 

In the following I will be primarily interested in the electromagnetic field, F^ v . Writ- 
ing this in spinor form, we have Faa'bb 1 with the antisymmetry ensuring that Faa'bb 1 = 
—Fbb'AA'- But using the above reduction, we note that this can be written as 

Fab = D Facbd'Za'b' + £ CD Fca'db'^ab) (18) 
I will use the notation that 
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Fab = -j^ C D Facbd' (19) 

so that 

Faa'bb' = Fab^a'B' + Fa'b'^ab (20) 
Maxwell's equations become 

V aa 'F A b = 4vr J B A ' (21) 

V AA 'F A , B , = AttJ a b , (22) 

where J^ 1 is the current source for the Maxwell field. The reality of and ensure that 

(Fab)* — Fa'b' and (Jab 1 )* = Jba 1 - 

Finally, if is another antisymmetric tensor, then 



F^G^ = 2(F AB G AB + F a ' b 'G A 'b>) (23) 
and if F and G are both real tensors, then 

F AB G ab = (F A ' B 'G A 'B'Y (24) 

Having established the notation, I will now state the theorem without proof [|T],0 . Given 
a spin s massless field without source, then the following integral gives the value of the field 
^(Q)ab...s at the apex (r=0) of the null cone emanating from a point in spacetime. 



T(0) AB - s mAB..s = / r S m(e)T AB ..si A i B ...i S o D o E ...o T (o x o Y 'V XY ^ DE .. T (25) 

Z7T J r,u const v 

-(2s + 1)^de..tO W l x o y 'Vxy'0 W S 



Here Tab..s is any covariantly constant spinor field, V xx'T A b.:s = 0. (In fact it need only 
be covariantly constant along the null cone of interest). This theorem states that the value 
of any massless field in flat spacetime can be determined by the integral over a sphere on 
the null cone emanating from that point. This expression is a generalisation to a dynamic 
massless field of arbitrary spin of the Kirkoff type integrals for static fields in terms of the 
integral over some surface surrounding the point in question of the normal derivatives of the 
field and the Green's function for that field. 

Although I have stated the theorem in terms of integrals over metric spheres on the null 
cone (r =const.), it can also be generalised to the integral over arbitrary two surfaces on the 
null cone. However, I will not use that generalisation here. 

While the use of this integral for the value of the field at points in the spacetime where the 
field is regular is interesting but unexceptional, a surprising result QXJ is that this integral 
also gives finite values if the point of interest is the location of a point charge (with its 
divergent Coulomb field). In fact, this integral (or rather the average of this integral over 
the future null cone emanating from the location of that charge at some time, and the past 
directed null cone emanating from that same point) gives exactly the radiation reaction field 
for an accelerating charge. 
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Frr = ~3^u (26) 

where the square brackets around the indices indicates anti-symmetrization, = ^(S^ u — 

S Vfl ). Ie, this Penrose integral automatically averages out the divergent field of the point 
particle to give just the finite radiation reaction contribution. Thus this Penrose integral 
approach differs substantially from techniques like the Dirac H or Abraham-Lorentz Q 
which suffer from divergences and the necessity for re-normalisations. 

The purpose of this paper will be to apply the above formula to the calculation of 
the radiation reaction field of an accelerating and time varying point dipole source for the 
electromagnetic field. Surprisingly, considering the fact the "Coulomb" portion of the field 
now diverges as 1/r 3 rather than as 1/r 2 for the case of a point charge, the Penrose integral is 
still finite, and gives a field which agrees with the radiation reaction field for a time varying 
unaccelerated dipole calculated by other methods |J. I will therefor assume that this finite 
field is also the correct radiation reaction field for an accelerated dipole. I have so far been 
unable to prove that this is consistent with the radiation damping one would calculate by 
more traditional techniques (eg, from the energy momentum tensor), but I can see no reason 
why it would not. 



II. SCALAR FIELD 

In order to gain practice, let me first calculate the radiation reaction field produced by 
a source which radiates scalar radiation. Consider a scalar field coupled to a point source 
travelling along the line x^(u) in flat spacetime, with a time varying source for the scalar 
field of intensity m(u), The Lagrangian for the massless scalar field is assumed to be 

S = J ^y^g<j) 4l <j) u d 3 x + Air J m(u)<j)(x(u))du (27) 

which gives the equation for <fi of 

V^V M = 4vr J m(u)5 4 (x fl - x{uY{u))du (28) 

In the urOip coordinates, the retarded Green's function for the scalar field in urOcf) coordinates 
is particularly simple, it is just 1/r. Thus the retarded solution for the scalar field is just 

0(n,r,^) = ^M (29) 
r 

We can now substitute this expression into the equation for the radiation reaction field. In 
order to write this in a slightly more transparent form, recall that l AA = o A o A . Thus 

o x o Y 'V<p= F<p^ (30) 
o a o x o y 'V X y'0 A = o A o A >o x o Y 'V X Y'{o A i A ') 

= Z M m a V Q m^ (31) 

since oao a = 0. This gives us, 
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,/ n 1 f / r\\ frh(u) .1 „. . .„... m(u). lm(u)\ , , 
0(0)= _/ r sin(0) - (- + r/(«) cos(0))(--^) - --^J rfdcty 

= -2m(u) (32) 

An interesting application is where the source for the scalar field is an internal oscillator 
with configuration variable q, such that m = eq. Ie, the Lagrangian is 

J {{- 9M,,9^x + J (±(q% - fiV ) + atf(s(u))) du (33) 

The equation of motion for the oscillator, including the effect of the radiation reaction field, 
is now 

- q, u , u - tfq + e(0 o (a;(u)) - 2eq, u ) = (34) 

where O is the value of the background field at the location of the particle. Ie, the radiation 
reaction field acts as a simple damping term to the internal harmonic oscillator, with damp- 
ing coefficient 2e 2 . In another paper I will use this to investigate the emission of radiation 
from an accelerating detector in interaction with the scalar field. 



III. ELECTROMAGNETIC RADIATION REACTION FIELDS 

To find the radiation reaction fields for the electromagnetic field, we must first solve the 
equations for the electromagnetic fields from a point dipole. First define the vector tangent 
to the path of the particle parallel transported over the null cone. 

T^ = l" + -n» (35) 

In addition, define the vectors 

Z»= (-F + \n») cos(#) - -i= sin(0) (m" + m") (36) 
2 y2 

X^= ^(-P + ^)sin(6) + cos(6)(m» + cos(0) + ^|(™ M - ™ M ) sin(0) (37) 

Y^= ^(-P + sin(0) + cos(0)(m" + m")^ sin(0) - -^(™ M - ™ M ) cos(0) (38) 

which are all vectors which are parallel over the whole of the surface of the cone u constant 
(but are not parallel off that cone.) For future needs, let me define the basis vectors 
such that 

(39) 
(40) 
(41) 
(42) 









= x» 


e (2)M 




e (3)M 


= 
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Define the dipole moment 

V»{u) = d x {u)X» + d y {u)Y» + d z {u)Z» = d {i) e {i ^ (43) 
as the dipole moment vector, with dm) zero. Furthermore, define 

= -(T^V U - T U V) (44) 

Then the vector potential for the electric dipole moment T> is 

A»(u,r,6,(l>) = W V S^ (45) 



with electromagnetic field 



Also define the tensor 



Ffju/ = V U A^ - V^A U (46) 



I Ztt 



= a {ij) e^ e W (47) 
where anj\ is antisymmetric in ij. Then the Penrose integral equation becomes 

J r sm(6)A A BL A L B oco D (o x o Y> V XY > F CD (48) 

— 3Fcdo w l x o x Vxx'Ow^j dOdcj) + C omplexC on jugate 

= -— [ rsin{e)Aa V rh tl n lJ l p rrf{l T V T F pa - 3F pa m a rh p V pl a )d9d4> + CC 
An J 

The first factor of half arises from the averaging over the future and past null cones (the 
contribution from the past null cone being zero). 

This integral, though very messy, can be evaluated, and gives the radiation reaction field. 
After extensive calculation, aided in an essential way with the GRTensorll computer algebra 
system 0, the result for the E and B fields are 

< = -5%^ - ^i^/ 2 («) + »«/(«)/(«) («) 

„« _\( D Fw a (j) v _ 2 D F\ya(j) D F w^(k) \ (i)(j)(k) ( 50 n 
RR 3 V Du 2 {k) Du Du ) y ' 

where D ™ is the Fermi Walker derivative of the quantity along the path of the particle 

-g^- = T»V u S a (5» - T»T a ). (51) 

and 

a (i) = e fT v V u T» (52) 
are the components of the acceleration in the XYZ frame. 
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We can find the radiation reaction fields for the case for a magnetic dipole moment, A4 M 
by simply taking the dual of all of the above equations, for which we get 



B 



RR — 



(0 _ 



2D 3 FW M^ D FW M® 
3 Du 3 Du 




(53) 



■W _ 



1 (Dl 



(54) 



E 



RR — 



We note that if the acceleration is not equal to zero, and the dipole moment is pe- 
riodic in time, the acceleration increases the radiation reaction field, since D 3 FW j^ ~ 
— Omega 2 Df wM ' ; unlike for the scalar field where the radiation reaction field was indepen- 
dent of the acceleration. 

One of the most interesting features of the radiation reaction field is that an accelerated 
magnetic dipole will have an electric component to its radiation reaction field. Thus, if the 
point dipole also has a non-zero charge, then the accelerated dipole will exert a force on 
that accelerated charge in addition to the normal radiation reaction force from the charge 
itself. On the other hand, since the radiation reaction field for an accelerated charge has no 
magnetic component in the rest frame of that charge, the accelerated charge will not alter 
the motion of the magnetic dipole moment. 

If the magnetic dipole moment is proportional to the angular momentum, via the mag- 
netic moment, \i of the point charge distribution, then the equation of motion for the dipole 
moment under constant magnitude of acceleration is 



This solution is subject to runaway behaviour just as is the accelerated charge. In this 
case the runaway behaviour is in the direction in which the angular momentum points, 
the angular momentum itself of course being conserved by the above equation. Choosing 
coordinates 8<j) such that the B points in the direction of 9 = 0, and the magnetic moment 
vector points in the direction O(w), $(«), and choosing g(u) = 0, F(u) constant, the 
equation of motion for O, $ are 



Using the lowest order solution for $(w), namely $ u = —fiB , and neglecting all but the 
lowest order terms in the equation for 0, we get 




(55) 



$(«),„= -fJ.B - 





(58) 



which has solution 



cos(O) = -tank -^-((/i5 ) 2 + f 2 )fiB (t - t ) 




(59) 
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However, retaining the higher order terms leads to runaway solutions, where in particular 
$ diverges exponentially. Ie, as usual, the electromagnetic radiation reaction is useful only 
in providing the lowest order corrections to the solution, and cannot be taken seriously as a 
complete solution. 



ACKNOWLEDGEMENTS 



I would like to thank P. Chen for inviting me to the Quantum Aspects of Beam Physics 
Workshop and J.D. Jackson for raising the issue of the radiation for circular accelerating 
beams and thus reviving my interest in this problem. I would also like to thank the Canadian 
Institute for Advanced Research and NSERC for support during this research. 



10 



REFERENCES 



[1] W. Unruh, Proc Roy. Soc. A348 447 (1976) 

[2] R. Penrose, W. Rindler Spinors and Space-time Cambridge U Press (Cambridge) 
(1984) 

[3] P.A.M. Dirac, Ann. de Flnst. H. Poincare 9 13 (1939) 

[4] J.D. Jackson Classical Electrodynamics, 2nd Ed. John Wiley (New York) (1975) 

[5] E.T. Newmann, R. Penrose J. Math. Phys 3, 566 (1962) 

[6] Daniel R. Stump, Gerald L. Pollack , Am. J. Phys. 65 81, (1997) 

[7] Peter Musgrave, Dennis Pollney, Kayll Lake, GRTensorll The GRTensorll package for 
Maple and Mathematica is available for free from |http: / / astro. queensu.ca~grtensor[ 



11 



